## [1] "ATTRIBUTE" "WETLAND_TY" "ACRES" "Shape_Leng" "Shape_Area"
## # A tibble: 9 × 2
## hpu result
## <fct> <dbl>
## 1 Bh 0.290
## 2 Bhs 0.118
## 3 Bhsm 0.168
## 4 Bi 0.227
## 5 E 0.149
## 6 H 0.583
## 7 I 0.259
## 8 O 0.0601
## 9 T 0.127
## Warning: Groups with fewer than two data points have been dropped.
## Warning: Removed 11 rows containing non-finite values (`stat_ydensity()`).
## Warning: Groups with fewer than two data points have been dropped.
## Warning: Removed 11 rows containing missing values (`geom_point()`).
Is there a model?
##
## Call:
## lm(formula = WET ~ hpu, data = hb_pedons_ext)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.28180 -0.12423 -0.06280 0.06577 0.71024
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.28980 0.02539 11.413 < 2e-16 ***
## hpuBhs -0.17225 0.04077 -4.225 3.08e-05 ***
## hpuBhsm -0.12180 0.07571 -1.609 0.10858
## hpuBi -0.06304 0.03787 -1.665 0.09689 .
## hpuE -0.14123 0.04905 -2.879 0.00424 **
## hpuH 0.29320 0.21547 1.361 0.17450
## hpuI -0.03110 0.05777 -0.538 0.59076
## hpuO -0.22975 0.05647 -4.069 5.88e-05 ***
## hpuT -0.16275 0.03300 -4.931 1.28e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.214 on 339 degrees of freedom
## (11 observations deleted due to missingness)
## Multiple R-squared: 0.1141, Adjusted R-squared: 0.09317
## F-statistic: 5.456 on 8 and 339 DF, p-value: 1.741e-06
## Analysis of Variance Table
##
## Response: WET
## Df Sum Sq Mean Sq F value Pr(>F)
## hpu 8 1.9984 0.249801 5.4562 1.741e-06 ***
## Residuals 339 15.5203 0.045783
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Warning: Removed 5 rows containing non-finite values (`stat_ydensity()`).
## Warning: Removed 5 rows containing missing values (`geom_point()`).
## `geom_smooth()` using formula = 'y ~ x'
## Warning: Removed 30 rows containing non-finite values (`stat_smooth()`).
## Warning: Removed 30 rows containing missing values (`geom_point()`).
Is there a model in W3?
##
## Call:
## lm(formula = WET ~ hpu, data = w3_pedons_ext)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.15737 -0.08136 0.00100 0.05010 0.22100
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.23275 0.03892 5.981 2.59e-06 ***
## hpuBhs podzol -0.06037 0.05504 -1.097 0.2827
## hpuBimodal podzol -0.14325 0.08702 -1.646 0.1118
## hpuE podzol -0.08355 0.06275 -1.331 0.1946
## hpuLithic Histosol 0.03725 0.06740 0.553 0.5852
## hpuTypical podzol -0.11575 0.06275 -1.845 0.0765 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.1101 on 26 degrees of freedom
## (5 observations deleted due to missingness)
## Multiple R-squared: 0.2335, Adjusted R-squared: 0.08605
## F-statistic: 1.584 on 5 and 26 DF, p-value: 0.1997
## Analysis of Variance Table
##
## Response: WET
## Df Sum Sq Mean Sq F value Pr(>F)
## hpu 5 0.095937 0.019187 1.5837 0.1997
## Residuals 26 0.315003 0.012115
##
## Call:
## lm(formula = WET ~ water_table, data = w3_pedons_ext)
##
## Residuals:
## 3 6 11 17 19 28 35
## -0.20724 -0.11312 0.04029 0.06486 0.19226 -0.04112 0.06408
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.350684 0.158615 2.211 0.078 .
## water_table -0.002206 0.002690 -0.820 0.450
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.1445 on 5 degrees of freedom
## (30 observations deleted due to missingness)
## Multiple R-squared: 0.1185, Adjusted R-squared: -0.05776
## F-statistic: 0.6724 on 1 and 5 DF, p-value: 0.4495
## Analysis of Variance Table
##
## Response: WET
## Df Sum Sq Mean Sq F value Pr(>F)
## water_table 1 0.014031 0.014030 0.6724 0.4495
## Residuals 5 0.104337 0.020867
Watershed 3 wells
## Warning: package 'stringr' was built under R version 4.3.1
Other soil moisture data
## plot WET
## 1 IL1 0.111505258
## 2 IL2 0.017028480
## 3 IL3 0.016389744
## 4 IH1 0.014945904
## 5 IH2 0.300747536
## 6 IH3 0.009245352
## 7 E01 0.108972447
## 8 E02 0.043977207
## 9 E03 0.042729688
## 10 E07 0.010488627
## 11 E08 0.022509384
## 12 E09 0.148845048
## 13 E10 0.342572154
## 14 E12 0.022831787
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 7674 rows containing non-finite values (`stat_smooth()`).
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 19943 rows containing non-finite values (`stat_smooth()`).
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 32704 rows containing non-finite values (`stat_smooth()`).
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 15780 rows containing non-finite values (`stat_smooth()`).
## Warning: Removed 32704 rows containing missing values (`geom_point()`).
## Warning: Removed 19943 rows containing missing values (`geom_point()`).
## Warning: Removed 7674 rows containing missing values (`geom_point()`).
## Warning: Removed 15780 rows containing missing values (`geom_point()`).
Old Denitrification soil moisture & oxygen
## ID UPL WET
## 1 1 0.9856686 0.01433141
## 2 2 0.9671750 0.03282499
## 3 3 0.9192316 0.08076836
## Warning: Removed 1 rows containing missing values (`geom_point()`).
## Removed 1 rows containing missing values (`geom_point()`).
## Removed 1 rows containing missing values (`geom_point()`).
## Warning: Removed 26 rows containing missing values (`geom_point()`).
## Warning: Removed 1 rows containing missing values (`geom_point()`).
## Warning: Removed 90 rows containing missing values (`geom_point()`).
## Warning: Removed 1 rows containing missing values (`geom_point()`).
## Removed 1 rows containing missing values (`geom_point()`).
## Removed 1 rows containing missing values (`geom_point()`).
## Warning: Removed 209 rows containing missing values (`geom_point()`).
## Warning: Removed 74 rows containing missing values (`geom_point()`).
## Warning: Removed 4989 rows containing missing values (`geom_point()`).
## `geom_smooth()` using formula = 'y ~ x'
## `geom_smooth()` using formula = 'y ~ x'